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Abstract 

A study is presented of central exclusive production of T(n5) states, where the 
T{nS) resonances decay to the final state, using pp collision data recorded 

by the LHCb experiment. The cross-section is measured in the rapidity range 
2 < y(T) < 4.5 where the muons are reconstructed in the pseudorapidity range 
2 < < 4.5. The data sample corresponds to an integrated luminosity of 2.9fb“^ 

and was collected at centre-of-mass energies of 7 TeV and 8 TeV. The measured 
T(15) and T(2S') production cross-sections are 

a{pp ^ pT{lS)p) = 9.0 ± 2.1 ± 1.7 pb and 
a{pp ^ pT{2S)p) = 1.3 ± 0.8 ± 0.3 pb, 

where the first uncertainties are statistical and the second are systematic. The T(IS') 
cross-section is also measured as a function of rapidity and is found to be in good 
agreement with Standard Model predictions. An upper limit is set at 3.4 pb at the 
95% confidence level for the exclusive T(3S') production cross-section, including 
possible contamination from Xbi^P) T(3S')7 decays. 
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1 Introduction 


Central exclusive production (CEP) of T{nS) [n = 1, 2, 3) resonances in pp collisions is 
thought to occur by photoproduction through the exchange of a photon and a pomeron (a 
colour-singlet system) between two protons, as illustrated in Fig. Since the protons do 
not dissociate, typically only a small component of momentum transverse to the beam 
direction (px) is exchanged in the interaction. The photoproduction of T resonances at 
LHCb can be computed using perturbative quantum chromodynamics (QCD), given the 
high photon-proton centre-of-mass energy, W, and the cross-section depends on the square 
of the gluon parton-density function, g(x), where Bjorken-x is the fraction of the proton’s 
momentum carried by the gluon [^. Measurements of the production cross-sections for the 
T(nS) resonances in the forward region covered by the LHCb detector are sensitive to g(x) 
in the region of small x down to approximately 1.5 x 10“®, where the knowledge of g(x) is 
limited. Furthermore, predictions for the T(nS) cross-sections at leading order (LO) and 
next-to-leading order (NLO) in the strong-interaction coupling differ greatly for the values 
of W probed in T(nS') resonance production, and there are signihcant variations depending 
on the models used to describe the T wave function and the f-channel exchange [l|^. 

Quarkonia photoproduction has been studied in exclusive production at HERA [4]-[^ , 
the Tevatron [lo] and the LHC |11 13 . At LHCb, exclusive production is associated 


with the absence of signihcant detector activity apart from that associated with the 
exclusive candidate. The background from proton dissociation occurring outside the 
detector acceptance is characterised as having a value of T candidate pr which is larger 
than that for exclusive production. 

In this article, the exclusive production cross-section of T(nS') resonances is measured in 
the hnal state where both muons lie in the pseudorapidity (g) range 2 < g{p,^) <4.5 
and the T{nS) candidate is reconstructed in the rapidity (y) range 2 < y{T{nS)) < 4.5. 
The pp data correspond to an integrated luminosity of 0.9 fb“^ at a pp centre-of-mass 
energy of y/s = 7TeV and 2.0fb“^ at a/s = 8TeV. Given the limited statistical precision, 
the data sets are combined to measure the production cross-sections. The LHCb detector 
and the simulated event samples are outlined in Sect. i In Sect. selection criteria 
are discussed, which exploit the absence of detector activity other than that associated 
with the T {nS) candidate. The signal efficiency and the various sources of background 


p 



Figure 1: Leading Feynman diagram for photoproduction of T(nS') states, where the photon- 
pomeron interaction is indicated by the shaded grey circle. 
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are also described. In Sect. |^two fits are described, which allow the determination of 
the exclnsive signal yield: by htting the T{nS) invariant mass spectrnm in order to 
separate T resonances from dimuon continuum background; and by htting the T(nS') 
candidate distribution to distinguish exclusively produced T resonances from those 
originating in hard interactions. Systematic uncertainties are summarised in Sect. and 
the measurements of the cross-sections are discussed in Sect.|^ Finally, for the T(IS') the 
differential cross-section, as a function of T(IS') candidate rapidity, is presented. 


2 Detector and simulation 


The LHCb detector 14,15 is a single-arm forward spectrometer, designed for the study 
of particles containing & or c quarks. It is fully instrumented in the pseudorapidity 
range 2 < rj < 5 and has tracking capability in the backward direction, in the range 
—3.5 < rj < —1.5. The detector includes a high-precision tracking system consisting of a 
silicon-strip vertex detector (VETO) surrounding the pp interaction region, a large-area 
silicon-strip detector located upstream of a dipole magnet with a bending power of about 
4 Tm, and three stations of silicon-strip detectors and straw drift tubes placed downstream 
of the magnet. The tracking system provides a measurement of momentum, p, of charged 
particles with a relative uncertainty that varies from 0.5% at low momentum to 1.0% at 
200 GeV/c. Photons, electrons and hadrons are identihed by a calorimeter system consisting 
of scintillating-pad (SPD) and preshower detectors, an electromagnetic calorimeter and a 
hadronic calorimeter. Muons are identified by a system composed of alternating layers of 
iron and multiwire proportional chambers. 

The trigger consists of a hardware stage, based on information from the calorimeter 
and muon systems, followed by a software stage, which applies a full event reconstruction. 
The hardware trigger requires events to contain at least one muon with a px greater 
than 200MeV/c. Low-multiplicity events are selected by requiring that fewer than ten 
hits should be detected in the scintillating pad detector, positioned just upstream of the 
electromagnetic calorimeter. In the subsequent software trigger, both of the hnal-state 
muons are required to have px > 400MeV/c. 

The exclusive production of T(nS') resonances is simulated using the SuperChiC 
software package [^, which provides the four-momentum of a single, transversely po¬ 
larised T(nS') resonance in each event. The decay of the T(nS') candidate is described 
by EvtGen [It], in which hnal-state radiation is generated using Photos 18 . The 


interaction of the generated particles with the detector, and its response, are implemented 
using the Geant4 toolkit as described in Ref. [^. Samples, each containing one 
million events, are prepared for T(?7,S') resonances decaying to In the same way, nine 

background samples of a similar size are prepared containing events where an exclusively 
produced Xfco,i, 2 (l-P, 2P, 3P) meson decays to the T(nS ')7 hnal states with T —)■ p+/i“. 
Separate samples are prepared for every Xb{mP) -P- T{nS)x (m, n = 1, 2, 3; n < m) decay. 
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3 Candidate selection 

Selection criteria are applied offline to events that pass the trigger requirements, to select 
well-reconstructed T(nS') candidates and to ensure the absence of unrelated detector 
activity. The latter set of requirements favours events containing a single pp interaction 
per bunch crossing. 

The hnal-state tracks, which must be associated with hits in the muon chambers, are 
required to lie in the pseudorapidity range 2 < <4.5 and to be of good quality. In 

extracting the differential cross-section for the T(IS'), the following intervals in T(IS') 
rapidity are considered; 2<y<3, 3<y< 3.5 and 3.5 < y < 4.5. Dimuon candidates 
are selected if the invariant mass falls in the range between 9 GeV/c^ and 20 GeV/c^, and 
the candidate p\ is less than 2GeV^/c^. The latter requirement favours photoproduction 
candidates, which have a characteristically low-pT- Events are rejected if one or more tracks 
are reconstructed in the backward direction. In the forward region exactly two tracks, 
corresponding to the muon candidates, are required, and these must be reconstructed both 
in the VELO and in the downstream tracking detectors. 

The selection criteria affect not only the T(nS') candidate but also the level of activity 
in the rest of the event, specifically through the requirements that there should be exactly 
two forward tracks, no backward tracks and fewer than ten SPD hits. The event is excluded 
if more than one proton-proton interaction occurs, causing a larger number of additional 
SPD hits or extra tracks to be reconstructed. The probability for an exclusive T event not 
to feature additional activity from another pp interaction in the same beam crossing is 
determined as the fraction of events containing no activity, according to these criteria, in 
a randomly accepted, hence unbiased, sample. After subtracting the contribution of the 
dimuon candidate, an event may contain fewer than eight SPD hits and no reconstructed 
tracks in the backward direction or tracks in the forward region. The fraction of randomly 
triggered events passing these criteria, /si, is found to be (23.63 ± 0.04)% for the 7TeV 
data and (18.48 ± 0.02)% for the 8TeV data. The difference arises because of the different 
beam conditions in the two data-taking periods, leading to a different average number of 
proton-proton interactions per event. 

The reconstruction, trigger and offline selection efficiencies are determined using 
simulated samples, and the combined efficiency varies between 77% and 84%. Eor signal 
candidates that pass the trigger and reconstruction stages, the offline selection criteria are 
more than 99% efficient. 

4 Determining the exclusive yield 

Candidates reconstructed in the 7TeV and 8TeV data sets are combined in a single 
sample, and two unbinned, extended, maximum-likelihood hts are carried out. A first ht 
is performed to the dimuon invariant mass spectrum, between 9GeV/c^ and 20GeV/c^. 
The ht contains a non-resonant background component and three resonant components. 
The three resonant components each receive contributions from exclusive signal, inelastic 
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Figure 2: Invariant dimuon mass spectrum for 7TeV and 8TeV data in the rapidity range 
2 < y(T) < 4.5 (black points). The fit PDF is superimposed (solid blue line). The T(15, 25, 35) 
signal components, used to derive weights, are indicated with a long-dashed (red) line, and the 
non-resonant background is marked with a short-dashed (grey) line. 


Table 1: Results of the invariant mass fits, within each rapidity interval. 


Parameter 

2 < 1/ < 4.5 

2 <y <3 

3<y<3.5 

3.5 <y< 4.5 

T(15,25,35) yield 
T(15) fraction 

T(25) fraction 

T(15) mass (MeV/c^) 

382 ± 26 

0.71 ± 0.03 

0.18 ± 0.03 

9452.5 ± 3.3 

146 ± 16 
0.74 ± 0.05 
0.16 ± 0.04 
9453.2 ± 4.3 

133 ± 16 

0.72 ± 0.06 

0.15 ± 0.05 
9452.4 ± 5.6 

94 ± 14 
0.68 ± 0.07 

0.26 ± 0.06 

9452.0 ± 9.0 


background and Xb feed-down decays. These contributions are indistinguishable in 

the invariant mass distribution. 

The probability density function (PDF) used to model each T{nS) signal peak is a 
Gaussian function with modihed tails (a double-sided crystal ball function |^). The mass 
differences for the T(25)—T(15) and T(35)—T(15) resonances are taken from Ref. |22] . 
The ratios of the T(25) and T(35) resolutions with respect to the T(15) are hxed to 
the ratio of their masses with respect to the mass of the T(15), following the procedure 
used in previous T measurements using LHCb data j^. The parameters that govern the 
shapes of the tails are taken from simulation, as is the resolution of the T(15) resonance, 
which varies from 35MeV/c^ to 57MeV/c^ in the different rapidity ranges. The yields of 
the signal components are all free to vary independently. 

A background PDF accounting for the non-resonant background is modelled using an 
exponential shape where the slope and normalisation are allowed to vary. 

The data are fitted in the whole rapidity range and in bins of rapidity. The ht results 
are given in Table HI and the £t in the full rapidity range is shown between 9 GeV/c^ and 
12 GeV/c2 in Fig. 
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Table 2: Estimated yields of feed-down background from —?■ T(nS ')7 decays in each 

T(nS') sample, where the uncertainties are statistical only. 


Signal window 

T sample 

Estimated contamination yield 
Xbi^P) Xbi‘2P) Xbi^P) 

2 < p(T) < 4.5 

T(l^) 

63 ± 10 

14 ± 5 

3 ± 2 


T(2^) 

— 

43 ± 12 

5 ± 3 


T(3^) 

— 

— 

21 ± 21 

2 < p(T) < 3 

T(15) 

31 ± 8 

2 ± 2 

0 ± 2 

3 < p(T) < 3.5 

T(l^) 

22 ± 6 

10 ± 4 

0 ± 2 

3.5 < y{T) < 4.5 

T(l^) 

8 ± 4 

0 ± 2 

3 ± 2 


Two sources of background contribute to the fitted signal: feed-down from x;, —)■ T 7 
decays, and inelastic interactions that involve the undetected products of proton dissociation 
or additional gluon radiation. 

The feed-down background is estimated using a combination of data and simulation, 
considering XbijnP) —t T(nS ')7 decays. Events are considered in the data set if exactly 
one photon is found in addition to the T candidate. Regions in the T 7 invariant mass 
spectrum are dehned, corresponding to the 2P, 3P) states, and the number of Xb 

candidates, for each decay Xft(mP) —)■ T(nS ')7 is counted. An estimate of the total 
feed-down content of the T data sample from each Xb state is found using the expression: 


iV, 


feed-down, 


^7 ^ ^mass-range 


( 1 ) 


Here P is the purity of the T {nS) in the corresponding mass window with respect to the non¬ 
resonant background, determined by fitting the dimuon mass spectrum for events 

with exactly one reconstructed photon; is the efficiency for reconstructing the photon 
produced in each Xbi'bnP) decay, determined using simnlated exclusive Xb{^P) 'T{nS)'y 
decays; and e 


mass-range 


= 0.9 corrects for the fraction of signal T candidates which are 
expected to fall outside the mass window. There are too few T( 3 S ')7 candidates to estimate 
the purity precisely so it is assumed to be 100%. Because of limited mass resolution and 
small sample sizes the Xb spin states cannot be resolved, so equal contributions from the 
Xfei(mP) and Xb 2 {'bnP) states are assumed. The Xbo radiative decay rate is expected to be 
relatively suppressed and is therefore neglected 
given in Table 

Since the mass shapes for signal and background do not significantly depend on px over 
the pt range considered, the Px distribution of the T candidates is determined using the 
A fit is then performed to the Px distribution, shown in Fig. nsing 


22 . The feed-down background yields are 


sPlot technique 24 


candidates in the full rapidity range 2.0 < p(T) < 4.5, with fit components corresponding 
to the T signal, inelastic background and feed-down background. The fraction of exclusive 
signal calculated from this fit is assumed to be the same for each rapidity bin. 

The Px distribution for the exclusive signal is derived from the simulated sample. 
At HERA, the distribution of the exclusive charmonium signal as a function of the 
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candidate was well described by an exponential fnnction, exp(— [^|^. The 
distribntion provides discrimination among varions prodnction sources because it 
approximates the squared four-momentum transfer, |f|, which depends on the production 
mechanism. Following Ref. 12 , Regge phenomenology is used to extrapolate the slope 


measured by HERA up to LHC energies according to the expression 


h{W) = 60 + 4a' log , 


( 2 ) 


where a' describes the slope of the exchange Regge trajectory and the constant 60 is specihc 
for interactions at a given photon-proton centre-of-mass energy, Wq [^. The SuperChiC 
generator models the pomeron-photon exchange and performs this extrapolation. 
Since the only published measurement of 6 x has very low precision j^, the generator is 
tuned to reproduce the LHCb measurement of bjj^ = 5.7 ± 0.1 GeV“^c^ 12 in exclusive 
J/'ifj production. The input values to SuperChiC are 60 = 5.6 GeV“^c^, a' = 0.2 GeV“^c^ 
and Wo = 90 GeV. The fit PDF is obtained from the simulated samples with these inputs, 
and a kernel estimation is employed to derive a shape to ht to data 25 . It is assumed 


that the inelastic background component is distributed according to a single exponential 
function |^. The slope and yield of this background function are free to vary in the fit. 

The total contamination from y;fc(mP) —)■ T(nS ')7 decays is constrained to be the sum of 
the contributions in Table and enters the p\ ht by means of a Gaussian constraint. Given 
that no analysis of exclusive Xb production has been undertaken, and the consequent lack of 
knowledge of the exclusive purity of the very small sample of reconstructed Xb candidates, 
it is assumed that inelastic processes contribute half of this feed-down background and 
the same inelastic background PDF is employed as that used to model inelastic T(nS') 
production. For the exclusive component, the shape of the dimuon PDF depends on 
the meson source, and the PDFs for each source, determined from simulation, are 

combined according to their relative contributions to the total feed-down yield. 

The exclusive purity, P, is dehned as the ratio of the exclusive signal yield, A^exdusive, to 
the number of candidates remaining in the sample, Ae^dusive + Aineiastic, after subtraction 
of the feed-down yield. The fit to the distribution in the full rapidity range, shown in 
Fig-i gives 

V = XT—-= (54 ± 11 ) %, 


A 


exclusive 


A 


inelastic 


with the exponential slope of the inelastic background measured to be — 0.21 ±0.26 GeV“^c^. 
The results of fits to the p^ distribution in each rapidity interval are consistent with this 
value. In order to validate the ht procedure, a set of pseudoexperiments is generated 
using the parameters obtained from the ht to the data, and the same ht is applied to each 
pseudoexperiment. The uncertainty on the purity is underestimated by 15% in the ht and 
the statistical uncertainty quoted takes account of this. 
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Figure 3: Fit to the p\ distribution of the T candidates in the full rapidity range. 


5 Systematic uncertainties 

The relative systematic uncertainties for the T(1S', 2S', SS") cross-sections in the various 
rapidity ranges are summarised in Table 

Contributions to the systematic uncertainty arising from the ht are considered: 
the uncertainty in the signal distribution as modelled by the SuperChiC generator 
and the variation of the exclusive signal PDF expected in the various rapidity bins. 
The SuperChiC generator is tnned to reprodnce measnrements of exclnsive J/'ip meson 
production made by LHCb [^. As no sufficiently precise measurements of the p\ 
distribntion in exclusive T(?7,S') resonance production exist, an estimate is made following 
Ref. Q, where it is argued from Regge theory that the slope 6o of the proton should be 
reduced by Aa'\og{m-r{nS)l'^j/ip)- A simulated sample is generated accordingly and used 
to derive a signal p\ template. Changing 6o from 5.6 to 4.7 produces a relative decrease in 
the exclusive yields of 6%, and this change is taken as the systematic uncertainty. For 
the differential cross-section measurements, the dependence of the signal p\ shape on 
rapidity is studied by replacing the exclusive signal p\ PDF with those determined in 
the smaller rapidity ranges, and the largest change in pnrity is taken as the nncertainty. 
Combining the systematic uncertainties in quadrature yields a total uncertainty for the 
exclnsive pnrity, V, between 7.2% and 8.2%. In addition, the possibility for variation in 
the shape of the continnnm dimnon backgronnd in p\ as a fnnction of mass is considered. 
The determination of the exclusive purity, R, is repeated in the dimuon invariant mass 
range from 9 to 12 GeV/c^, and the difference is taken as a conservative estimate of the 
systematic uncertainty. In Table these sources contribute to the uncertainty denoted 
‘pnrity fit’. 

The uncertainty arising from the p\ shape derived from simulation and used to 
describe the feed-down background is considered separately. The feed-down backgronnd 
PDF is constrncted nsing only contribntions from the Xbi(mP) and Xb 2 {'mP) backgronnd 
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components, in equal parts, assuming no contribution from xwijnP) decays. Since it is 
not possible to resolve the spin states in data, we consider a conservative change where 
the nominal PDF is replaced with that for background originating from the decay of a 
Xboi'^P) meson. The ht to data is repeated and the change in exclusive purity is taken as 
the associated uncertainty. In addition, there is uncertainty associated with the exclusive 
fraction of the Xbi'^P) feed-down background, which in turn affects the overall shape of 
the feed-down background used in the spectrum. Since the size of the data set is too 
small to allow a data-driven estimate of the Xb(mP) sample exclusive purity, the PDF for 
this purity is assumed to be uniform between 0 % and 100 %, and the effect of changing it 
by ±1 standard deviation (0.50 ± 0.29) is therefore considered. The resulting change in 
the exclusive T(nS') sample purity is taken as the systematic uncertainty. The estimate of 
the yield of XbijnP) meson feed-down in data includes the determination of the photon 
reconstruction efficiency using simulated samples of Xfe(mP) —)■ T(? 7 ,S ')7 decays. Samples 
of —)■ and —)■ decays are used to validate the agreement 

between photon reconstruction efficiencies in data and simulation. An uncertainty of 
5% is taken as the systematic uncertainty on the photon reconstruction efficiency to 
account for the small differences seen [^. The resulting uncertainty is very small for the 
T(IS') but larger for the T(2S', 3S') samples where the relative contamination from the 
X;,(mP) background is larger. These three systematic uncertainties on the cross-section are 
combined in quadrature and are presented as the ‘feed-down b.g.’ systematic uncertainty 
in Table [H 

An estimate of the contamination of the T(IS') and T(2S') samples from the decays 
T(25) —>■ T(lS') 7 r° 7 r°, T(3S') —>■ T( 2 S'){ 7 r° 7 r°, 77 }, is made using the observed T(nS') 
candidate yields in data and the relevant T' —)■ TX and —)■ fp'fiT branching fractions 

1^ . The estimated contaminations are taken as systematic uncertainties for the T(IS') 
and T(2S') cross-sections. 

Uncertainties in the PDFs used to fit the T(nS') candidate invariant mass spectrum are 
considered. Alternative T(nS') signal PDFs are produced, obtained using kernel estimation. 
The systematic uncertainties on the exclusive purity and each of the yields are assessed 
using pseudoexperiments generated with the nominal invariant mass PDFs and htted with 
a model where the signal PDFs are replaced by those obtained using kernel estimation. 
The effect of replacing the exponential PDF used to model the non-resonant background 
in the invariant mass fit with a second-order polynomial function is found in the same way 
using pseudoexperiments. Combining these two sources of uncertainty in quadrature leads 
to a relative uncertainty which is less than 4% for all the cross-sections. This uncertainty 
is labelled ‘mass fit’ in Table [H 

The LHCb integrated luminosity has been measured with a relative uncertainty of 1.7% 
at 7TeV and 1.2% at 8 TeV 27 . The integrated luminosity is multiplied by the estimated. 


selection-dependent, fraction of events that contain no interactions other than the one 
that produces the signal candidate, /si. The determination of /si from data depends 
upon the subtraction of the T signal candidate’s SPD hits. The spread in the signal 
candidate SPD hit multiplicity is estimated from data to be one hit, and the fraction /si 
is therefore recomputed with the signal subtraction increased from two to three and the 
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Table 3: Summary of the relative systematic uncertainties, in %. 



2 <y <3 
T{1S) 

3<y<3.5 

T{1S) 

3.5 <y < 4.5 
T(1R) 

2 < 1/ < 4.5 

T(1R) T{2S) T{3S) 

Purity fit 

14.2 

14.2 

14.2 

13.7 

13.7 

13.7 

Feed-down b.g. 

12.2 

12.2 

12.3 

12.2 

14.6 

12.5 

T' feed-down 

4.0 

4.3 

5.4 

4.5 

11.1 

— 

Mass fit 

2.2 

2.8 

2.9 

2.1 

2.8 

3.6 

Luminosity 

2.3 

2.3 

2.3 

2.3 

2.3 

2.3 

B{T^p+p-) 

2.0 

2.0 

2.0 

2.0 

8.8 

9.6 

Total 

19.5 

19.7 

20.0 

19.3 

24.8 

21.4 


change is taken to be the systematic uncertainty. To account for variations as a function of 
data-taking time, the variation of the estimated single-interaction fraction is evaluated in 
each uninterrupted period of data-taking during which conditions are constant, typically 
an hour long, instead of considering each year as a whole, and the change with respect to 
the nominal fraction is taken as the systematic uncertainty. Combining the uncertainties 
in quadrature yields an overall relative uncertainty for each year of 2.3%. The systematic 
uncertainties on the luminosity for each year are assumed to be 100% correlated. 

The branching fractions, i3, for T(IS'), T(2S') and T(3S') to decay to the dimuon 
final state are accounted for to determine the T(nS') production cross-section. These 
branching fractions are taken from Ref. 22 and, for the T(IS'), T(2S') and T(3S') states. 


carry relative uncertainties of 2%, 9% and 10%, respectively. These are propagated to an 
uncertainty on the production cross-section. 


6 Cross-section 

The cross-section is obtained using 


N, 


a = 


exclusive 


C X e X B(T{nS) —)■ /i+/i“)' 


(3) 


The effective integrated luminosity, C, is 580 pb taking into account the values of /si 
given in Sect. [3 27 


The quantity e is the efficiency correction, which is obtained in each rapidity bin and 
for each resonance, and which is averaged for 7TeV and 8TeV data according to the 
luminosity in each year, and B{T{nS) —)■ is the T(?7,S') —>■ branching fraction. 

The measured exclusive production cross-sections in the LHCb acceptance are 

cr(pp —)■ pT(lS')p) = 9.0 ± 2.1 ± 1.7 pb, 

a{pp pT{2S)p) = 1.3 ± 0.8 ± 0.3 pb, and 

a{pp — pT{3S)p) < 3.4 pb at the 95% confidence level. 
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Table 4: Production cross-section for the T(IS') resonance in ranges of T(IS') rapidity, where the 
muons are required to lie in the pseudorapidity range 2 < < 4.5. The first uncertainties 

are statistical and the second are systematic. 


2<y<3 3 < 2 / < 3.5 3.5 < r/ < 4.5 

ct(T( 15)) (pb) 3.4 ± 0.9 ± 0.7 2.9 ± 0.8 ± 0.6 2.6 ± 0.8 ± 0.5 


Table 5: Measured da{T{lS))/dy, where the data have been corrected for the effect of the LHCb 
geometrical acceptance. The statistical and systematic uncertainties are combined in quadrature. 

2<y<3 3 < 2 / < 3.5 3.5 < r/ < 4.5 
da{T{lS))/dy (pb) 8.8 ± 3.0 7.8 ± 2.7 7.1 ± 2.6 


where the hrst uncertainties are statistical and the second are systematic, and where 
the limit on cr(T(3S')) includes possible contamination from Xb feed-down. The limit is 
calculated using pseudo-experiments and includes the effect of systematic uncertainties, 
where correlations are assumed to be negligible. The T(IS') production cross-section is 
given in smaller ranges of T(IS') rapidity in Table 

After correction for the LHCb geometrical acceptance, the cross-sections in Table can 
be compared to theoretical predictions. The efficiency for an T candidate to be produced 
in the range 2 < y{T) < 4.5, and to decay to muons which lie inside the acceptance, 
2 < < 4.5, is 45%. In the smaller ranges of T(IS') rapidity considered for the 

differential cross-section measurement, the efficiency is lowest in the outer ranges, at 39% 
{2 < y < 3) and 36% (3.5 < y < 4.5), and highest in the central range, at around 74% 
{3 < y < 3.5). The correction depends on the rapidity distribution in the simulated sample, 
which has a different shape to those of, for example, the predictions in Fig. |^. To estimate 
the systematic uncertainty on the geometrical acceptance correction, the simulated samples 
are reweighted to obtain a uniform rapidity distribution within each rapidity bin, and 
the change in the geometrical acceptance is taken as the systematic uncertainty. This 
corresponds to a relative change in the geometrical acceptance of less than 6%. The 
differential cross-sections, da{T(IS))/dy are given in Table and are shown in Fig. 
compared to LO and NLO predictions Q. The LHCb data are in good agreement with 
the NLO prediction. 

The LHCb results are also compared to theoretical predictions for the underlying 
photon-proton cross-section, as a function of the centre-of-mass energy of the photon- 
proton system, IT, as shown in Fig. |4 )d. There are two contributions to the photoproduction 
of an T(nS) resonance, depending on which proton emits the virtual photon. The pp 
cross-section is given by 

d.-(pp^pT(lg)p) ^ ^ 

where the predictions for the photon-proton cross-section are weighted by absorptive correc¬ 
tions S‘^(W±) and the photon fluxes ^ for photons of energy k± (Mx(nS)/2) exp(±||/|). 
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Figure 4: Measurements of exclusive T(IS') photoproduction compared to theoretical predictions. 
In (a), the T(IS') cross-section in bins of rapidity is shown, compared to LO and NLO predictions. 
The LHCb measurements are indicated by black points with error bars for uncorrelated errors, 
and solid rectangles indicating the total uncertainty. In (b), the photon-proton cross-sections 
extracted from the LHCb results are indicated by black points, where the statistical and 
systematic uncertainties are combined in quadrature. The entire kF-region in which these LHCb 
measurements are sensitive is indicated. Measurements made by HI and ZEUS in the low-kF 
region are indicated by red and blue markers, respectively [^[^[^. Predictions from Ref. are 
included, resulting from LO and NLO fits to exclusive J/'0 production data. The filled bands 
indicate the theoretical uncertainties on the 7 TeV prediction and the solid lines indicate the 
central values of the predictions for 8 TeV. In (b) predictions from Ref. using different models 
for the T(15) wave function are included, indicated by ‘bCGC’. 


The absorptive corrections and photon fluxes are computed following Ref. [^. 

The three bins of T(IS') rapidity chosen in this analysis correspond to ranges of W 
for the LU+ and LU_ solutions. The contribution to the total cross-section from the IU_ 
solutions is expected to be small and is therefore neglected. The dominant W+ solutions 
are therefore estimated assuming that they dominate the cross-section, and are shown in 
Fig. |4j). The magnitude of the theoretical prediction for the LF_ solutions is added as 
a systematic uncertainty. The good agreement with the NLO prediction seen in Fig. 
is reproduced. The LHCb measurements probe a new kinematic region complementary 
to that studied at HERA |^[^[^, as seen in Fig. |^, and discriminate between LO and 
NLO predictions. In Fig. |^, the LHCb data are also compared to the predictions given in 
Ref. using models conforming to the colour glass condensate (CGC) formalism that 


take Into account the f-dependence of the differential cross-section. All agree well with the 
data. The solid (black) and dotted (blue) lines correspond to two different models for the 
scalar part of the vector-meson wave function. 
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7 Conclusion 


The first measurement of exclusive T{nS) production in pp collisions at 7 and 8TeV is 
presented, and a differential cross-section is extracted as a function of T(IS') candidate 
rapidity. The data probe a previously unexplored kinematic region in photon-proton centre- 
of-mass energy. The results are compared to theoretical predictions and a strong preference 
for those including next-to-leading order calculations is seen. Exclusive production studies 
at LHCb will be improved during LHC Run II following the installation of scintillators at 
high IpI, which will allow for improved trigger efficiency for exclusive production processes 
and additional suppression of the background from inelastic interactions 
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